Preparation and detection of d-wave superfluidity in two-dimensional optical superlattices 
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We propose a controlled method to create and detect d-wave superfluidity with ultracold fermionic atoms 
loaded in two-dimensional optical superlattices. Our scheme consists in preparing an array of nearest-neighbor 
coupled square plaquettes or "supeiplaquettes" and using them as building blocks to construct a d-wave super- 
fluid state. We describe how to use the coherent dynamical evolution in such a system to experimentally probe 
the pairing mechanism. We also derive the zero temperature phase diagram of the fermions in a checkerboard 
lattice (many weakly coupled plaquettes) and show that by tuning the inter-plaquette tunneling spin-dependently 
or varying the filling factor one can drive the system into a d-wave superfluid phase or a Cooper pair density 
wave phase. We discuss the use of noise correlation measurements to experimentally probe these phases. 



Ultracold atoms in optical lattices|[T]| are promising simu- 
lators of complex many-body problems and model Hamilto- 
nians that arise in condensed matter physics. They provide a 
clean system with parameters which can be tuned in a con- 
trolled fashion from the weakly interacting to the strongly in- 
teracting limits. The observation of the superfluid to Mott in- 
sulator transition with bosons|2 1 and recent experimental real- 
ization of both repulsive and attractive Hubbard models||3] ID 
with Fermions are important steps in that direction. Recently, 
the existence of superexchange antifeiTomagnetic correlations 
have been shown in an array of isolated double wells|5j [SJ. 
Here we extend similar considerations for a d-wave superfluid 
state and propose an experimental scheme to realize and de- 
tect it. 

The repulsive Hubbard model on a square lattice is one of 
the most important model Hamiltonians for strongly interact- 
ing fermions and is widely believed to contain the essential 
features of high temperature superconductors|7|. At half fill- 
ing (one fermion per site), strong repulsion leads to localiza- 
tion of the fermions resulting in a Mott insulator. Superex- 
change interactions lead to an antiferromagnetic ground state. 
Away from half filling, antiferromagnetic coiTelations com- 
pete with the kinetic energy of the holes and the phase diagram 
remains unknown. Although there is no rigorous proof of ex- 
istence of superfluidity in this model, approximate theoretical 
and Monte Carlo simulations at high temperatures [8J have 
conjectured the presence of d-wave superfluid ground states 
away from half-filling. 

It has been known for some time|9| that the minimal 2D 
unit that can sustain d-wave pairing physics is a single square 
plaquette (see Figjl]). The single plaquette has been the start- 
ing point of various approachesi 10 11 12, 13| to the Hubbard 
model on the square lattice. Even with purely repulsive inter- 
actions, there is a range of parameters where two holes tend 
to bind together on a single plaquette rather than to delocalize 
among different plaquettes. The hole pair that is created has 
a d-wave symmetry and can lead to d-wave superfluidity once 
the plaquettes are coupled. 

We first propose an experimental scheme to verify these 
concepts by creating an array of isolated super-plaquettes (i.e. 




C) I I I I I I I 




FIG. 1: A plaquette is the minimum system that exhibits d-wave 
symmetry, a) When loaded with four fermions the ground state is d 
wave symmetric while when loaded with 2 the ground state exhibits 
s wave symmetry. Consequently the two states have non zero matrix 
element with the d-wave pair creation operator. Here we consider the 
situations when plaquettes are coupled into a super-plaquette array 
and into a checkerboard array, which are schematically represented 
in b) and c) respectively. In the picture the spin independent intra- 
plaquette tunneling J is represented by a thick solid line and the 
inter-plaquette tunneling, ^ J by a dashed line. The subscript a 
in emphasizes that it can depend on the spin of the atom. 



two adjacent plaquettes coupled by a weak tunneling, see Fig. 
[TJi) and loading them with six Fermions each. To measure 
the binding energy of the hole pair, we analyze the coher- 
ent dynamics of the super-plaquette. Finally we derive the 
phase diagram of the weakly coupled plaquettes with spin- 
dependent inter-plaquette tunneling. We show that the sys- 
tem can be driven to either a d-wave superfluid or a density 
wave state |fT3). We propose an experimental scheme to realize 
this and discuss the detection of the two quantum phases via 
atomic noise correlations measurements lfT4l [TSl [161 [TTl [TSl . 
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FIG. 2: a) The red solid line corresponds to the pair binding energy in a plaquette. For < U / J < 4.6, Af, > and consequently it is 
energetically favorable to have two holes in the same plaquette. b) Coupling parameters as function of [//J for the effective XXZ Hamiltonian, 
Eq. ipl. The blue line corresponds to g^/A;, as a function of U /J. This is the only parameter that appears in the super-plaquette Hamiltonian. 
When more than two plaquettes are coupled, new virtual processes have to be accounted for which change the Ising term coupling constant. 
While the latter become proportional to (?f /At, the transverse coupling constant remains the same, i.e. proportional to g^/Aj,. The red line 
shows the absolute value of the difference between these two parameters. 



I. PLAQUETTE FERMION MODELS 



Consider fermions in an isolated plaquette (shown in 
Fig|T^). Assuming one accessible single particle state in each 
well (i.e. level spacing much larger than other energy scales in 
the problem), the system is described by the Hubbard Hamil- 
tonian 

H = -J ^ cl^c^,„ + U^n^rni,r (1) 

{r,r'),cr r 

where J is the tunneling matrix element and U is the onsite 
Hubbard repulsion. Here c^^ are fermionic annihilation oper- 
ators, hra = cj^^c^g. are number operators, r = 1, ... 4, labels 
the four sites in a plaquette and the term (r, r') indicates that 
the sum is restricted to nearest neighbors. 

The eigenstates in a single plaquette depend on the filling 
factor. When filled with = 4 or = 2 fermions, the 
ground state is a spin-singlet exhibiting d and s wave symme- 
try respectively. We denote the corresponding ground states 
as |4) and |2) (See Fig[l}i). On the other hand, for iV = 3, 
the ground state is degenerate with = 1 /2 and ± ipy 
symmetry in the regime U < Ut ^ 18.6 J. We denote them 
as IS^"^'^-*) with a and t specifying the spin polarization and 
the orbital "chirality" (r ~ ±). The d\s s wave symmetry of 
the |4) and |2) states is the crucial element in obtaining the d- 
wave pairing mechanism, and the hole-pair creation operator 
that connects the two states must have a d-wave symmetry |9l; 



i.e. (2|A^|4) ^ 0, where 

Ad = (Si2 + S34 - Sl4 - S23)/2, (2) 

and Srr' ~ {cl^cl,^ — cl^cl,^)/\/2 creates a singlet on the rr' 
bond. 

We next look at the case of 2 holes in two isolated plaque- 
ttes. The holes can bind together within the same plaquette 
or separate as single holes in each plaquette depending on the 
binding energy defined as 

Ai, = 2Eg{N = 3) - Eg{N = 4) - Eg{N = 2) (3) 

being positive or negative respectively. Here Eg{N — n) is 
the single plaquette ground state energy when loaded with n 
atoms. As shown in Fig. |2j Af, is a non-monotonic function 
ofU/J ifTTl IT2I [131 [T9I , which reaches a maximum value of 
Af, K, 0.04J at f/ « 2.45J and becomes negative for U = 
Uc > 4.58 J. Consequently only when U < Uc hole pairs on 
a single plaquette are energetically stable. 

Two adjacent plaquettes can be coupled through a weak 
(possibly spin dependent) tunneling to form a super- 
plaquette (see Fig{T]). As long as < ^ A;,, the states 
|4, 2) and |2,4) are lower in energy and occupation of the 
states are energetically suppressed. They can 

only be populated as "virtual" intermediate states leading to an 
effective super-exchange interaction between |4,2) and |2,4). 
Specifically, by treating the |4) and |2) states as the pseudo- 
spin components I fl") and I -IJ-) of an effective spin 1/2 system, 
the interaction between the effective spins can be described by 
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FIG. 3: Loading of fermions in a patterned loading array of plaquettes: a) First a unit filled band insulator of fermionic atoms is prepared in 
a rectangular lattice with lattice period 4/3A and 2 A along the x and y direction respectively, b) A phase shifted 4A lattice along the x axis 
is slowly ramped up. c)As the period 2 A lattice is ramped up, the period 4/3 A is ramped down generating a double well lattice. The bias is 
chosen such that the ground state has four fermions in the right and two fermions in the left sites of each double well, d) The lattice depth of 
the 2A lattice is increased to make all wells independent. Subsequently the bias is turned off and the isolated wells split into two. The overall 
scheme creates an array of independent plaquettes loaded in the ground state with alternating filling factors of 2 and 4 along the x direction. 



an XXZ-type Hamiltonian 



H, 



2Ah 



(4) 

with (J°L^'^'^ standard Pauli matrices acting on the right (R) 
or left (L) effective pseudo-spins and g the coupling matrix 
element between the right and left plaquettes which is of order 
one (see Fig|2]). 

Eq.(|4|i contains the essential physics we are interested in 
this work. If J| — J| the energy eigenstates are effective 
triplet and singlet states \t,s) = ^(|4,2) ± |2,4)) with the 
triplet being the ground state. These are separated by an en- 
ergy gap ^ Ab from the rest of the Hilbert space. The ground 
state has a non zero expectation value of the d-wave pair cor- 
relation operator (t|A^AL|t) ^ which leads to a d-wave 
superfluid when coupling all the plaquettes. On the contrary 
if J| ^ J|, the Ising term dominates and any infinitesimal 
symmetry-breaking perturbation will collapse the state into 
|4,2) or |2,4) which are inherently density-ordered states. 
These considerations indicate that the many body phase dia- 
gram of this model would depend on the spin-dependent cou- 
plings in a non-trivial way. 



II. PREPARATION AND DETECTION 

To verify the energy structure of the hole-pair states we 
want to create an array of isolated super-plaquettes loaded 
with 2 and 4 fermions in the left and the right plaquettes re- 
spectively. 

An aiTay of plaquettes can be created by superimposing two 
orthogonal optical super-lattices formed by the superposition 
of two independent sinusoidal potentials which differ in pe- 
riodicity by a factor of two |6, 20 21 1. The aim here is first 
to load the fermions in a 2D array of independent plaquettes 
with alternating filling factors of 4 and 2 along one direction. 



The preparation procedure we propose is based on a unity- 
filled band insulator |,22J. and uses adiabatic manipulations of 
a super-lattice potential created by standing waves with four 
different periodicities A, 2 A, 4A and 4A/3 (For other strate- 
gies see lE l). Such wavelength combinations are available 
for typical Fermionic atoms such as ^"K and ^Li or can be 
engineered by intersecting four pairs of laser beams with ap- 
propriate angles |23|, resulting in a set of four equidistant k- 
vectors. For the discussion we will also assume that there is 
a deep axial lattice that freezes the atom motion along the z 
direction. 

The process of patterned loading an array of plaquettes is 
depicted in Figj3] The initial band insulator is formed in a 
rectangular lattice with period 4/3 A and 2 A along the x and 
y directions, respectively. Adiabatically introducing a phase 
shifted 4A lattice along the x axis creates an effective 4A- 
lattice with three sites per cell with different energy offsets 
(See Figj3|l. By ramping down the period 4/3A lattice while 
ramping up the period 2A one can convert the potential into a 
double well lattice with an asymmetric energy offset. This 
offset has to be large compared to the on-site repulsion to 
guarantee that the final ground state corresponds to a system 
with four fermions in the lower and two fermions in the higher 
sites of each double well. After increasing the lattice depth 
of the 2A lattice to suppress tunneling between the wells the 
bias can be removed and the isolated wells can be split into 
four sites by slowly turning on the A lattice along both x and 
y directions. As a result, an array of independent plaquettes 
in the ground state with alternating filling factors of 4 and 2 
along the x direction is created. To obtain an array of super- 
plaquettes, the tunneling between adjacent plaquettes can be 
controlled by the depth of the 2A potential, while the 4A lat- 
tice isolates super-plaquettes. The extra 4A/3 lattice is needed 
to balance the offset created when the long lattice is added. 
Spin dependent control of the inter-plaquette tunneling can 
be achieved by additional control over the laser polarizations 
I121J. 

Under coherent quantum evolution, this system will exhibit 
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FIG. 4: Time evolution of the population of super-plaquettes in 
Ni2{t) (blue), N2A{t) (purple) and N33{t) (yellow) configurations 
respectively. The fast frequency component is determined by Aj, and 
the slow envelope by J'^y^/Afc. Here U/J = 2.6 and J'/ J = 0.01. 

Rabi oscillations between |4, 2) and |2, 4) states with the |3, 3) 
populated virtually (as seen in Fig. Hll. The frequency of oscil- 
lation of the envelope is given by 4^ J'^/Af, and can be used 
to measure the binding energy of the hole pairs. These ideas 
are similar to the ones used to measure the super-exchange 
energy in an array of double wells|5|. We next describe a 
method to measure the number of super-plaquettes in |4,2), 
1 2, 4) and |3, 3) states using band-mapping techniques. 

The distribution of super-plaquettes in different states can 
be obtained from the number distribution of fermions on the 
left and right plaquettes separately. To obtain the number dis- 
tribution, after quenching J', each plaquette is first converted 
to a single site by removing the short lattice potential. By 
counting the number of the atoms in each Brillouin zone af- 
ter a band-mapping sequence on release from the lattice, the 
average occupation of all bands is determined. In a second 
step, the sequence is repeated, but the vibrational level of each 
atom on the right plaquette is increased by one after the con- 
version to a single well. This can be achieved in analogy to the 
scheme demonstrated in |24|, by splitting the combined well 
and recombining in the presence of an energy bias. The bias is 
created using a phase shift on the 4-lambda lattice in such a 
way that only one of the plaquettes is affected. By comparing 
the resulting band occupations with those of step 1, the aver- 
age filling of each band on each side is reconstructed. Both 
of these steps can be repeated with an additional Feshbach 
resonance sweep which converts all atom pairs to molecules. 
The remaining single atoms allow the separation of the frac- 
tions of even and the odd number states in each band of each 
side, which completes the determination of the plaquette atom 
number statistics. 

From the obtained histogram of atom number distributions 
for the left and right plaquette at time t, Ni2{t) , N24,{t) and 



NQ3{t) can be determined ll25l to measure the hole-pair bind- 
ing energy. 

III. PHASE DIAGRAM ON CHECKERBOARD LATTICE 

Now we study the more general case of a checkerboard ar- 
ray of plaquettes such as the one shown in FigjT];. We will 
restrict our analysis to the regime |Ab| ^ gj^ and Af, > 
where we can treat the states |2) and |4) as the low energy 
modes and adiabatically eliminate high energy states via sec- 
ond order perturbation theory. This procedure yields a more 
general effective XXZ Hamiltonian given by 

Heff = E [--^^(^R^R' + ^R^R') + ^'^R'^R'] (5) 
<R,R'> 

where = ^^3^. = i^'i^ + Jf)^,- Here we have 
neglected terms proportional to <t^ since it is a conserved 
quantity in these systems. The coupling matrix element 
has a complicated dependence on U/ J as shown in Figj2] For 
details about its derivation see methods. 

The zero temperature phase diagram of the XXZ Hamilto- 
nian is known and consequently can be used to infer the phase 
diagram of the corresponding fermionic system lfT3l l26l l27l . 
At 3/8 fermionic filling, i.e. iV4 = N2, there is a second 
order phase transition as the J^- / ratio is varied: While 
for J^- 1 < 1 the ground state corresponds to a gapless 
d-wave superfluid ( a magnetically ordered phase in the XY- 
plane for the effective spins), for J^- / > 1 it becomes a 
gapped Cooper pair density wave state, d-pdw (an antiferro- 
magnetic order phase for the effective spins). The pdw state 
is not the usual particle-hole charge density wave state, but 
can be viewed as a crystal of d-wave Cooperpairs |28|. The 
point = J-^ is the critical point. In Fig. pa, we show the 
phase diagram as a function of J| / J| where the tendency of 
anisotropic tunneling to stabilize the d-pdw phase can be ob- 
served. For the spin independent case J| = J|, our phase 
diagram is in agreement with the one obtained in Ref. lfTSl . 
exhibiting a critical point at Ucr ^ 2.7 J. Notice that since 
our analysis is based on the assumption that gj' ^ A;,, the 
parameter regime where it is applicable considerably reduces 
as one approaches the points U = and U/J^ 4.6 where 
A;, vanishes. 

Away from 3/8 fermionic filling, the phase diagram is al- 
most insensitive to the J| / J| ratio. There is a first order phase 
transition from the d-pdw phase to the d-wave superfluid as 
the chemical potential is varied away from 3/8 filling. Due 
to the first order character of the transition, the d-pdw phase 
is surrounded by a small region where one observes phase co- 
existence 11131 127I . Except this region and the state at special 
filling factors: 1/4, 3/4, 1/2 and 1 when the system turns into 
an insulator, the low energy phase is always a gapless d-wave 
superfluid as shown in Fig.jSj). 

To explore the different quantum phases, we propose to first 
prepare an initial array of isolated plaquettes, with the desired 
filling factor (see methods), and then slowly increase the inter- 
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FIG. 5: a) Zero temperature phase diagram at 3/8 fermionic filling (A^4 — N2 — Ntotai/'^)- At this filling the system exhibits a second order 
phase transition from a d-wave Superfluid, d-SF, (magnetic ordered state in the XY plane for the effective spins) to a d-pair density wave state, 
d-pdw, (antiferromagnetic ordered phase for the effective spins ) when the axial, J^, and the transverse J^, coupling constants that appear in 
the XXZ Hamiltonian become equal. As both and J' depend on the inter-plaquette tunnelings J| | and the ratio U / J, the critical point can 
be controlled by tuning these microscopic parameters, b) Zero temperature phase diagram of the effective XXZ Hamiltonian as a function of 
the filling factor: Away from 3/8 fermionic filling = N2), the d-wave superfluid phase is energetically favorable except at the fermionic 
filling factors 1/4 (or = Ntotai ) and 1/2 or (A'4 = Ntotai ) where the system becomes an insulator. There is an additional first order phase 
transition from a d-pdw to a d-Sf as the filling factor is varied away A'4 = and this induces phase separation in a small range of filling 
factors around A'4 = N2). 



plaquette coupling taking advantage of the experimental abil- 
ity to tune the lattice geometry and the ratio J^/ J-'- by using 
spin dependent control. The effective XXZ Hamiltonian al- 
lows us to estimate the entropy of the ordered state (near Tc) 
and hence the T/Tp of the initial trapped Fermi gas without 
the optical lattice required for observing these phases (assum- 
ing all later modifications are adiabatic). This gives an esti- 
mate of T/Tf ~ 0.01 to observe the d-wave superfluid and 
T/Tf ^ 0.03 to observe the charge ordered state. These val- 
ues should be achievable with current technology. 



IV. DETECTION OF THE QUANTUM PHASES VIA NOISE 
CORRELATIONS 

Noise correlations lfl4l can be used to detect the different 
quantum phases. The atomic noise in an expanding cloud is 
related to the following four point functions at the time of the 
release: 

Gq5j, cx (nQo-nQv) - {n(ia){nQ'a') (6) 

with fiQa cx X); s ('^L'^scr), being the quasi- 

momentum distribution and Lj^ a vector connecting the lattice 
sites s and I. G'^'^ (Q', Q) has the required ingredients to dis- 
tinguish the two phases. On one hand, G^^(Q,Q') contains 
terms proportional to |(A)p (see Eq|2|, and consequently a 
d-wave superfluid with (A) ^ must exhibit interference 
fringes at Q + Q' = Kn, K being the reciprocal lattice vector 
of the plaquette array, which is half of the reciprocal lattice 



vector of the underlying lattice. The d wave nature of the 
state will be signaled by a modulation of the peaks with an a 
overall envelope with the characteristic d-wave nodal planes 
along Qx — ±Qy and Q'^ = ^Q'y shown in Fig|6] In- 
formation about the density order is given by G^^ which will 
show sharp dips at Q — Q' — {2n + l)K/2 in the presence 
of a d-pdw phase and a flat profile in the superfluid phase. In 
contrast to the unequal spin correlations, G^^ will not have 
the inherent d — wave symmetry. Equal spin correlations will 
also always exhibit dips at 5{Q — Q' — nK) which reflect the 
characteristic anti-bunching of fermions |17|. Unfortunately, 
the amplitude of these dips can be about 30 times stronger 
than the dips at (2?i + l)K/2 inherent from the charge density 
order (see Fig. This caveat can be avoided by quenching 
the kinetic energy of the atoms in the plaquettes before the re- 
lease by slowly merging each plaquette into a single well. For 
details about the noise interference pattern we refer the reader 
to the method section. 



V. CONCLUSIONS 

In summary we have described a technique to prepare and 
detect d-wave superfluidity in ultra-cold fermionic systems 
loaded in optical super-lattices. Even though our theoretical 
predictions are restricted to a regime where perturbative ex- 
pansions are valid, the ideas presented in this letter might open 
a window to explore in the laboratory more complex regimes 
which do not allow for a clear theoretical description. We 
point out that so far we have ignored the presence of any any 
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FIG. 6: The left panels display the unequal spin noise correlations, G^^{Q, Q' = Q + mk), where k = 2Tv/ax is the reciprocal lattice vector 
of the underlying lattice with spacing a, assuming a d— wave superfluid. The right panel shows equal spin noise correlations (with opposite 
sign) — G^^ (Q, Q' = — Q) for a d-pdw state. Here we have subtracted the local density-density correlations within a plaquette (see methods). 
G^^(Q, Q') exhibits interference peaks at Q + Q' = Km in the d-wave superfluid regime with K = k/2 is the reciprocal lattice vector of the 
plaquette array. The peaks are modulated by an overall envelope which probes the d wave symmetry. This modulation causes the disappearance 
of the peaks along the nodal lines Qx = ^Qy On the contrary — G^^ (Q, — Q) exhibits interference peaks at Q — Q' = K(2m + l)/2 in the 
d-pdw state. Due to the large kinetic energy of the fermions within a plaquette these peaks however are much weaker than the anti-bunching 
peaks, which always appear at Q + Q' = nK. 



inhomogeneous confining potential in the checkerboard lattice 
case. As the binding energy is a very sensitive quantity to any 
bias within a plaquette, its presence can certainly prevent the 
observation of the predicted phases. Consequently in this type 
of experiments it would be ideal to use a "box" type trapping 



potential such as the one proposed in Ref. llSTI 
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Methods 

Beyond the Effective Hamiltonian in a plaquette 

The validity of the effective Hamiltonian, Eq.( |4]i can 
be checked in the super-plaquette system by going one 
step beyond and deriving a more general effective Hamil- 
tonian within the Hilbert subspace spanned by all the di- 
rect products of the low energy eigenstates of the isolated 
plaquette i.e |4,2), |2,4) and eight possible IS^*^'^', 3^'*'^') 
configurations. Even though in total there are 10 states 
only 4 of them are coupled by J^. These are \t), 
1+) = E..(-l)'^^-^|3('^'^\3(<^.^))/2, \s) and h) = 
J2aT |3('^'^\ 3('^'^))/2. In this basis the generalized effective 
Hamiltonian reduces to 



H = 








-gJ' 































-qSJ' 


V 








-gdJ' 


A, 



(7) 



Here 5 J' = j'T-j'T and J' = j'T + j'i. The Hamiltonian is 
block diagonal so the eigenstates are independent linear com- 
binations of {\t), 1+)} and {|s), |— )} respectively and have 



energies: 



2,4 



(8) 



For gJ'^ ^ |Ab| the eigenstates are grouped in two dou- 
blets with common intra-doublet splitting Ag^J ^ J /Ai, and 
doublet splitting ~ A^. The sign of A;, determines which 
one of the two doublets is higher in energy. Since first or- 
der tunneling transfers amplitudes across the two doublets, if 
gJ'^ ^ \Ab\, these processes are energetically costly and only 
second order processes within the doublets are relevant. This 
assumption yields the effective spin Hamiltonian, Eq.(|4]i. On 
the contrary for gj'^ ^ |Ab| there is no clear energy separa- 
tion between the doublets and the ground state is a mixture of 
the various independent plaquette states. 

The corrections from the effective Hamiltonian can be mea- 
sured by probing the coherent dynamical evolution proposed 
in the main text. In the weak coupling regime N2i{t) and 
N42{t) will exhibit nice sinusoidal oscillations with a fre- 
quency proportional to g^J'^/Ai,. However for A^ ^ J' a 
more complex dynamics will occur. More quantitatively, us- 
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ing Eq.(|7]l one gets that 



1 



2 sin^[a;i3i/2] 



13 



± 



2uj 



13 



(9) 



'13 



Here = uji ~ u)^. In the Umit Ab 3> J', dominates 
over LOi, the |3, 3) states are only virtually populated and the 
above expressions reduce to 



N2i,i2{t)^ ^{^±cos{iu^t)) N33{t)^0 (10) 



Derivation of tlie Effective Hamiltonian parameters: 
Multi-plaquette case 

The microscopic parameters of the many-plaquette Hamil- 
tonian are more complicated that the ones in Eq.(|4|i since 
in addition to the virtual couplings between |4i^,2i{+i) <^ 
\'^rAr+i) thi-ough tunneling to |3fl,3ij+i) (which are in- 
versely proportional to Af,), when many plaquettes are cou- 
pled, there are additional processes that one must consider: 
We have to include the couplings between |4fl,4fl+i) <^ 
\^Ri^R+i) mediated by |5ij,3fl+i) states and between 
|2fl,27?,+i) ^ \2r,2b.+i) mediated by \3r,1r+i) (which are 
inversely proportional to 174 /A44 and 5I/A44 where A44 = 
E5 + E3 — 2E4, A22 = E3 + El — 2E2 and 54^2 are cor- 
responding coupling matrix elements). Since the energy gaps 
A44 and A22 are about one order of magnitude larger than Af,, 
to be consistent, one has in addition to sum over all possible 
higher energy virtually populated states (leading to additional 
corrections we denote by 5). All these processes give rise to 
an effective coupling constant along the axial direction given 



by: gl 
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2A22 



2A44 



Initial preparation of tlie plaquettes with N2 — 2N4 

The preparation starts by loading a unit filled band insulator 
of fermionic atoms in a rectangular lattice with lattice period 
3/2 A and A along the x and y direction respectively. After 
this, the A lattice is ramped up along x simultaneously with a 
weaker 3A lattice. While the former creates an additional well 
in between two initially populated wells, the latter provides an 
energy offset to guarantee that the new wells remain unpopu- 
lated. Next, one has to slowly turn off the 3/2 A lattice while 
turning on a 2A lattice. At the end of this procedure one has 
generated a distribution of double wells along x with the ap- 
propriated bias such that the ground state population has 242 
pattern. After increasing the lattice depth of the 2A lattice un- 
til all wells are made independent one can suddenly turn off 
the 3A lattice which provided the bias without causing any ex- 
citation. The overall scheme creates an array of independent 



plaquettes loaded in the ground state with alternating filling 
factors of 2, 4, 2 along the x direction. 



Noise Correlations 

Here we provide explicit expressions for the noise correla- 
tions shown in Fig|6] These expressions can be derived using 
the fact that the states to probe are linear superpositions of 
states populated with only 4 and 2 fermions in each plaquette, 
and therefore that they can be mapped to an effective spin 1 /2 
system. Using this simplification, the noise correlations can 
be written as 

GTi(Q,Q') (X A(Q,Q') + S+-(Q + Q')I?(Q,Q') dD 
G"(Q,Q') « A(Q,Q')-F+(Q,Q')<5(Q-Q'-toK) 
-5^^(Q-Q')F-(Q,Q') (12) 



The term A corresponds to trivial diagonal density-density 
correlations within a plaquette and ideally one could subtract 
this term from the measured noise pattern. (Q + Q') = 



E 



,i(Q+Q')-(R- 



is the analog of the quasi- 



momentum distribution when one maps the spins to hard- 
core bosons and consequently probes the long range order 
If the system is a d — wave superfluid (Q + Q') will 
exhibit sharp interference peaks at Q + Q' = Kn, K be- 
ing the reciprocal lattice vector of the plaquette array, which 
is half of the reciprocal lattice vector of the underlying lat- 
tice. On the contrary 5*^ (Q + Q') will have a flat pro- 
file in the case of a d-cdw state. The term I'(Q,Q') — 



pi(Q-l,.+Q'l,') 



(2|4^c:^|4) 



probes non-diagonal local 



plaquette correlations (r, r' label sites within the same pla- 
quette) and gives information about the d— wave nature of 
the superfluid. X'(Q,Q') modulates the S+-(Q + Q') peak 
structure inducing the characteristic d-wave nodal lines along 
Qx = ±Qy and Q'^ = iQy, where the signal vanishes due to 
the d-wave symmetry of the Cooper pairs (see Fig|6]l. 

Information about the density order is given by G^^ . 
S'^^(Q - Q') = J2R^R^e''^^-^'^<^-^'^al^al^,) shows sharp 
peaks at Q = Q' + {2n + l)K/2 in the presence of a d-cdw 
phase and a flat profile in the superfluid phase. In contrast 
to the unequal spin correlations the terms -F*(Q,Q') do not 
have the inherent d — wave symmetry. They are just geomet- 
ric factors that depend on the kinetic energy within the isolated 
plaquettes: i^±(Q,Q') = |a4(Q, Q')±-a2(Q, Q')IV4 with 

«n(Q,Q') = Err' e^^^''""*^' ^ , and n - 2,4. 

The delta function 6{Q — Q' — mK) reflects the characteristic 
anti-bunching of fermions. 



